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ABSTRACT 

We have modelled the inner surface brightness profiles of 39 alleged 'core' galaxies with the core- 
Sersic model, and provide new physical parameters for the largest ever sample of 'core' galaxies fit 
with this model. When present, additional nuclear components were simultaneously modelled and the 
typical rms scatter of the fits (out to ~10") is 0.02 mag arcsec -2 . Model-independent estimates of 
each core's break radius are shown to agree with those from the core-Sersic model, and a comparison 
with the Nuker model is provided. We found an absence of cores in what amounts to 18% of the 
sample which are reclassified here as Sersic galaxies with low values of n (< 4) and thus shallow inner 
profile slopes. In general, galaxies with n < 3 and a < 183 km s _1 do not have depleted cores. We 
derive updated relations between core-Sersic break radii, their associated surface brightness, bulge 
luminosity, central velocity dispersion, and predicted black hole mass for galaxies with depleted cores. 
With the possible exception of NGC 584, we confirm that the inner negative logarithmic profile slopes 
7 are < 0.3 for the 'core' galaxies, and > 7 > —0.1 for six of these. Finally, the central stellar mass 
deficits are found to have values typically within a factor of 4 of the expected central black hole mass. 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: fundamental parameter — galax- 
ies: nuclei — galaxies: photometry — galaxies:structure 



1. INTRODUCTION 

The stellar distributions in galaxies have played a valu- 
able role in guiding our understanding of the galax- 
ies themselves. In particular, the accessibility of high- 
resolution imaging offered by the Hubble Space Telescope 
(HST) substantially advanced our appreciation of the 
complexity of galaxy cores (e.g. Crane et al. 1993; Kor- 
mendy et al. 1994; Jaffe et al. 1994; Ferrarese et al. 1994; 
Grillmair et al. 1994; van den Bosch et al. 1994; Lauer 
et al. 1995; Byun et al. 1996; Gebhardt et al. 1996; Car- 
olio et al. 1997; Faber et al. 1997). For instance, the 
centers of real galaxies may contain such distinct com- 
ponents as bright active galactic nuclei (AGN), nuclear 
star clusters, flattened nuclear discs and bars, dust lanes 
and clouds. On the other hand, giant stellar evacuation 
zones are also observed. Luminous galaxies with such 
shallow cores had of course long been known to exist 
from ground-based observations (e.g. King & Minkowski 
1966, 1972; King 1978; Young et al. 1978; Binney & Ma- 
mon 1982; see the review by Graham 2012a) but HST 
enabled us to accurately quantify these. 

After studying 14 bright elliptical galaxies with the 
pre-refurbished ffST/WFPCl, Ferrarese et al. (1994) in- 
troduced a 4-parameter double power law model to de- 
scribe the inner surface brightness distributions of bright 
galaxies. While the (relatively brighter) galaxies in their 
sample which possessed depleted cores with shallow inner 
profiles were grouped as "Type I" , the remaining galax- 
ies, labeled "Type II" , had a profile that remained steep 
all the way into the center. Examining a larger sample 
of galaxies imaged using the same HST /WFPC1 high- 
resolution Planetary Camera, Kormcndy et al. (1994) 
and Lauer et al. (1995) largely agreed with the division of 
galaxies presented in Ferrarese et al. (1994), referring to 
them as 'core' galaxies and 'power-law' galaxies, respec- 
tively. They also advanced a double power law model 



which they dubbed the 'Nuker law' for fitting the (un- 
derlying host galaxy) surface brightness profiles of early- 
type galaxies. The Nuker model had an additional fifth 
parameter to moderate the transition between the two 
power laws — as introduced by Hcrnquist (1990, his eq. 
43) for modelling the internal density profiles of galaxies. 

The physical process (es) responsible for the observed 
difference between the inner surface brightness profiles 
of 'core' galaxies and the fainter 'power-law' galaxies 
(nowadays referred to as 'Sersic' galaxies because these 
spheroids have Sersic light profiles rather than power-law 
light profiles) provide valuable clues about the galaxies' 
past history. In bright galaxies, the widely advocated 
'dry' (i.e. gas poor) galaxy merger hypothesis (e.g. Faber 
et al. 1997) can result in the gravitational sling shot of 
central stars (core scouring) due to the coalescence of 
supermassive black holes (SMBHs) from the progenitor 
galaxies (e.g. Begelman, Blandford & Rces 1980; Makino 
& Ebisuzaki 1996; Merritt & Milosavljevic 2005; Merritt 
2006). It is possible that the sizes and mass deficits of 
such partially depleted cores may reflect the amount of 
merging and damage caused by the black holes (after 
having eroded any pre-existing nuclear stellar compo- 
nents: Bekki & Graham 2010). Having an accurate quan- 
tification of the physical parameters defining the centers 
of galaxies is therefore important. Moreover, reliable 
break radii Rb, used to denote the sizes of the cores, may 
even be useful for predicting black hole masses (Lauer et 
al. 2007a). 

While investigating the lack of any connection be- 
tween the double power-law model and the curved galaxy 
brightness profiles observed outside of the cores, Graham 
et al. (2003, see their figures 2-4) revealed that the Nuker 
model's break radius, and other parameters, were not ro- 
bust quantities but are sensitive to the radial range of 
the surface brightness profile that is fitted. For example, 
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the break radii Rb were shown to vary by more than a 
factor of three. The parameters' sensitivity was recog- 
nised to arise from the Nuker model's efforts to fit an 
outer power-law to what is actually a curved brightness 
profile. The luminosity profiles of bright (core) galaxies 
(Mb < —20.5 mag) — which show a downward devia- 
tion from the inward extrapolation of their outer Scrsic 
(1963, 1968) profile — were subsequently shown to be 
precisely represented by the core-Sersic model (Graham 
et al. 2003; Trujillo et al. 2004). 

While Lauer et al. (2005) missed this development, 
Ferrarese et al. (2006) found the core-Sersic model to 
be highly applicable to bright early-type galaxies in the 
Virgo cluster. Lauer et al. (2007a, b) subsequently wrote 
that "Graham et al. (2003) have criticized the Nuker rb 
as being sensitive to the domain over which the Nuker 
law was fitted, particularly when the outer limit of the 
fit extends only slightly beyond r&. In practice, however, 
the Nuker laws are fitted over a large radial range that 
extends well beyond r&." However this was not the prob- 
lem identified by Graham et al. (2003), who had demon- 
strated that the Nuker model parameters deviated fur- 
ther from the true values as the fitted radial extent was 
increased. 

Based on "work in preparation" Lauer et al. (2007b) 
refuted that their Nuker break radii were biased "in any 
way" because their radii reportedly agreed very well with 
model-independent values of where the curvature in the 
surface brightness profile was a maximum. This was a 
surprising claim because these latter values should not 
be dependent on the radial extent of the data while the 
Nuker model break radii are a strong function of the fit- 
ted radial extent (Graham et al. 2003). Kormendy et al. 
(2009, their section 4.1) subsequently buoyed the Nuker 
model and dismissed the core-Sersic model. Giiltckin et 
al. (2009) then over-looked any and all concerns about 
the Nuker model which they presented along with Nuker 
model parameters from Lauer et al. (2005), and encour- 
aged readers to use this data, additionally noting that the 
surface brightness profiles that were fit with the Nuker 
model are available at the Nuker web page 1 . Giiltckin et 
al. (2011) continued in this vein, motivating us to further 
investigate, nearly a decade on, the issue of whether the 
Nuker model parameters are reliable, physically mean- 
ingful quantities, or if instead the core-Sersic model pa- 
rameters may be preferable. At stake is not only the 
accuracy to which we quantify the cores of galaxies, but 
our subsequent understanding of cores and how they re- 
late to their galaxy at large. 

In this paper, we focus on the nuclear structure of 
galaxies by re-analyzing the surface brightness profiles 
of all 39 'core' galaxies imaged with the WFC2 / F555W 
or F606W filter and listed in Lauer et al. (2005) to be 
a 'core' galaxy (see section [2]). For reference, Trujillo et 
al. (2004) modelled only 9 possible 'core' galaxies, Fer- 
rarese et al. (2006) modelled 10, and Richings, Uttley 
& Kroding (2011) have very recently modelled 21 'core' 
galaxies. We are therefore modelling the largest sample 
of suspected 'core' galaxies to date. For comparison's 
sake with the Nuker model break radii, we use exactly 
the same surface brightness profiles as Lauer et al. (2005), 
available at the previously mentioned Nuker web-page. 

1 http:/ /www. noao.edu/noao/stafT/laucr/wfpc2_profs/ 



We first concentrate on measuring the core size using 
the core-Sersic model (see sections [3] and @| . We ad- 
ditionally take the nuclear excess, usually nuclear star 
clusters or AGN emission, into account while modeling 
the underlying host galaxy light. In section [5] we use two 
model-independent core size estimators and reveal that 
one of these can not be used while the other is consistent 
with our core-Sersic break radii. Furthermore, we con- 
firm that the published Nuker model break radii are typi- 
cally 100% larger than the break in the surface brightness 
profile determined relative to the inward extrapolation of 
the outer Scrsic function (Trujillo et al. 2004). We ad- 
ditionally report that 'artificial' break radii have been 
reported in what were alleged to be 'core' galaxies but 
are actually Sersic galaxies with no break in their Scrsic 
profile and which thus have no partially depleted core 
relative to their outer light profile (section |6]) . Through- 
out this paper we use terms such as 'actual', 'true' and 
'real' break radii and cores when referring to galaxies 
that have inner surface brightness profiles which break 
downward from (i.e. have lower flux than) the inward 
extrapolation of the outer Sersic model which describes 
their outer stellar distribution. 

Sets of structural parameter relations encompassing 
central as well as global properties are presented in sec- 
tion [71 In particular, equations involving the break ra- 
dius and associated surface brightness, and the lumi- 
nosity, are derived. We also investigate the (core size)- 
(central black hole mass) relation in section 17.21 Using 
updated data, we find that the break radius can be used 
to consistently predict the black hole mass when using 
either the Mbh~<J or Mbh~L relations for 'core' galaxies 
when coupled with our updated Rb-o and Rb-L relations. 
We go on to discuss the detection of additional nuclear 
components in the full sample of Scrsic and core-Sersic 
galaxies in section [8] while section [9] summarises our main 
conclusions. 

2. SAMPLE SELECTION 

Lauer et al. (2005) analyzed and presented fits to the 
major-axis surface brightness profiles of 77 relatively 
bright, nearby, early-type galaxies. Every galaxy in 
their sample was observed with the Wide Field Plane- 
tary Camera 2 (WFPC2: Biretta et al. 2001) onboard 
the HST and was centered on the PC CCD (which has 
an image scale of 0".0456 pixel^ 1 and a field view of 
800 x 800 pixels). Their sample lacks any characterizing 
selection criteria and comprises galaxies at distances of 
~ 10-100 Mpc. While almost all of the 'core' galaxies 
were imaged with the F555W filter (similar to broad- 
band V), two were not and we use their F606W data 
(roughly broadband R) instead. We refer the reader to 
Lauer et al. (2005) for an extensive description of the 
sample and images, which includes procedures adopted 
for PSF deconvolution 2 , dust obscuration correction, sky 
subtraction and background source masking. 



Our sample comprises 39 early-type galaxies which are 
the 'core' galaxy subset of the 77 galaxies presented in 

2 Some of the merits and disadvantages of image deconvolution 
are described in Ferrarese et al. (2006). 
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Fig. 1. — Left panel: A 3-parameter Sersic plus 2-parameter point-source model describes the entire available radial extent of NGC 596 
while the 5-parameter Nuker model (fit taken from Lauer et al. 2005) fails to describe the inner and outer light profile. Right panel: A 
3-parameter Sersic plus 2-parameter point-source model describes NGC 1426 better than the 5-parameter Nuker fit taken from Lauer et 
al. (2005) which can not describe the nucleus. The rms scatter, A, about the major-axis, F555W-band surface brightness profiles pertains 
only to the data points included in the fit (shown by the filled symbols). 



Lauer et al. (2005). This sample selection enables a di- 
rect comparison with the (published) Nuker model's es- 
timation of the core size and related parameters, and to 
further achieve this direct comparison we have modelled 
the same published light profiles 1 . The global properties 
of our target galaxies are summarized in Table 1, which 
presents their morphology, magnitude, updated distance 
and velocity dispersion. For the six lenticular galaxies 
plus one Sa spiral galaxy, we have roughly converted 
their total galaxy magnitudes, reported by Lauer et al. 
(2007b), into bulge magnitudes using a mean T^-band 
bulge-to-disc ratio of 1/3, equivalent to a mean bulge-to- 
total ratio of 1/4 (Graham & Worley 2008; Laurikainen 
et al. 2010, their section 6.3 and references therein). 
The la range on this B/T ratio for SO-Sa galaxies is 
about a factor of 2, corresponding to a la uncertainty 
of ±0.75 mag for our bulge magnitudes. In passing we 
note that for four of these seven disk galaxies (NGC 507, 
NGC 2300, NGC 3607, and the Sa galaxy NGC 7213), 
the dust-corrected X-band B/T ratios reported by Lau- 
rikainen et al. (2010) are 0.33, 0.28, 0.33 and 0.18, re- 
spectively, suggesting that we are not too far off with 
our adopted i?-band value of 0.25. Laurikainen et al. do 
however identify NGC 3706 as an elliptical galaxy, while 
de Vaucouleurs et al. (1991) refer to it as SAO(rs), and 
we have no comparison B/T ratios for NGC 4382 and 



NGC 6849. In any event, it should be remembered that 
these latter three galaxies make up a fairly small fraction 
of the total sample. 



3. THE SERSIC AND CORE-SERSIC MODELS 

By combining CCD images of early-type Virgo cluster 
galaxies with deep, large field-of-view, photographic im- 
ages, Caon et al. (1993) revealed that the Sersic (1963) 
model fits the main parts of the profiles of both el- 
liptical and spheroidal galaxies astonishingly well over 
large ranges in surface brightness (down to ^29 B-m&g 
arcsec -2 ). In general, after excluding the nuclear region, 
the Sersic R 1 /™ model fits the surface brightness profile 
of both elliptical galaxies and the bulges of disk galaxies 
remarkably well over their entire radial range (Caon, Ca- 
paccioli & D'Onofrio 1993; D'Onofrio, Capaccioli & Caon 
1994; Young & Currie 1994; Andedakis et al. 1995; Gra- 
ham et al. 1996). In fact, systematic deviations from this 
model appear to signal cither the presence of a central 
light deficit or an additional nuclear component (Balcells 
et al. 2003; Graham & Guzman 2003). Indeed, using 
CCD images with the extended profiles from Caon et al. 
(1993) and other photographic data, one of main conclu- 
sions noted in Kormendy et al. (2009) was exactly this. 

The radial intensity distribution of the 3-paramctcr 



4 



Dullo & Graham 



TABLE 1 

Updated global parameters of the 'core' galaxy 
sample from lauer et al. (2005). 
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Notes. — Col. (1) Galaxy name. Col. (2) Morpho- 
logical classification from the NASA/IPAC Extragalactic 
Database (NED)". Col. (3) Absolute K-band (galaxy or 
bulge) magnitude obtained from Lauer et al. (2007b) and 
adjusted using the distance from col. (4). Sources: (t) 
Tonry et al. (2001) after reducing their distance moduli 
by 0.06 mag (Blakeslee et al. 2002); (n) from NED (3K 
CMB). Col. (5) Central velocity dispersion from Hyper- 
Leda' (Paturel et al. 2003). The superscript * is used to 
indicate two 'power-law' galaxies taken from Lauer et al. 
(2005) and used for illustrative purpose in Fig. HI 

"(http:/ /nedwww. ipac.caltcch.edu ) 

f> (http:/ /leda. univ-lyonl.fr) 

Sersic R}l n model, a generalization of de Vaucoulers 
(1948) 2-parameter i? 1 / 4 model, is defined as 



I(R) = I e exp < -b n 



(1) 



where I e denotes the intensity at the half-light radius R e . 
The quantity b n ~ 2n — |, for 1 < n < 10 (e.g. Caon et 
al. 1993), is a function of the shape parameter n, and is 
defined in a way to ensure that R e encloses half of the 
total luminosity. A review of the Sersic model and its 
associated expressions can be found in Graham & Driver 



(2005). 

Many, but not all, spheroids fainter than -20.5 B-mag 
contain additional nuclear components. This is illus- 
trated in Fig. Q] where we display a Nuker model and 
a Sersic plus point source model fit to the surface bright- 
ness profiles of NGC 596 and NGC 1426. The pro- 
files have been taken from the Nuker web-pages, and 
the fitted Nuker model parameters are from Lauer et 
al. (2005). Both models agree on the absence of a par- 
tially depleted core, and as such the galaxies are classi- 
fied as Sersic galaxies. There are, however, considerable 
differences regarding the quality of the fits. For NGC 
596, the 3-paramctcr Sersic plus 2-parameter Gaussian 
model accommodates the entire observed radial range 
of the brightness profile remarkably well with a smaller 
root mean square (rms) residual than the 5-paramctcr 
Nuker model. The Nuker model fit (taken from Lauer et 
al. 2005) to this galaxy's light profile not only excluded 
the inner most data points but also clearly reveals a sig- 
nificant departure from the profile in the outer region. 
Similarly, while the Sersic model plus Gaussian function 
can represent the entire observed profile of NGC 1426, 
the 5-parameter Nuker model cannot describe the extra 
compact light source at the center and has more scatter 
in the residual profile. 

As mentioned previously, the light profiles of luminous 
(Mb < —20.5 mag) elliptical galaxies depart systemati- 
cally from the Sersic model near their center. It is impor- 
tant to realize that this departure, a downward deviation 
with respect to the inward extrapolation of the outer 
Sersic profile, emanates from a central starlight deficit 
and is not due to dust (which would result in a dramatic 
color change). Such stellar distributions can be described 
using the core-Sersic model introduced by Graham et al. 
(2003) and applied in Trujillo et al. (2004). A blend of 
an inner power-law and an outer Sersic function, it can 
be written as 



7/0 



exp 



jja _|_ jja \ !/H 



R? 



with 



/' = / fe 2-^ a exp \b(2 1 / a R b /R e ) 1 / n 



(2) 
(3) 



lb is intensity at the core's break radius Rb, 7 is the 
slope of the inner power law region, and a controls the 
sharpness of the transition between the inner power-law 
and the outer Sersic profile. As in the Sersic model, 
R e is the effective half-light radius of the outer Sersic 
function and b has the same general definition as before 3 . 
In practice, the 6-parameter core-Sersic model can be 
reduced to a 5-parameter model by setting a to some 
large constant value. Trujillo et al. (2004) set a — > 00, 
so that the transition from Sersic profile to power law 
at Rb is infinitely sharp, with no transition region, an 
approach effectively adopted by Ferrarese et al. (2006). 
In this paper, as in Richings et al. (2011), we explore and 
use a range of finite values for a. 

4. FITTING ANALYSIS 

3 One recovers the Sersic R 1 /™ function from equation (2) when 
setting Rj, and 7 to zero. 
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We fit the one dimensional light distributions of the 
39 underlying host galaxies using Sersic and core-Scrsic 
models, and account for any additional nuclear light com- 
ponents with either a Gaussian or an exponential model. 
Fits to the full 39 galaxies are available in Appendix A. 
The quality of the fits, as indicated by the rms given in 
each panel of Appendix A Fig. 21, is excellent, typically 
0.01-0.03 mag arcsec~ 2 . 

As can be seen, the 3-parametcr Sersic model prof- 
fers a good match to two (NGC 4473 and NGC 5576) of 
the 39 galaxies 4 all the way to the HST resolution limit. 
While these galaxies have shallow central profiles, the 
profiles do not 'break' from the outer envelope — fur- 
ther evidenced by the small values of a used with the 
Nukcr model (Lauer et al. 2005). Application of the 
Sersic model, along with small inner Gaussian and ex- 
ponential functions, yields a satisfactory fit to the lu- 
minosity profiles of an additional five nucleated galaxies 
(NGC 1374, NGC 4458, NGC 4478, NGC 4486B and 
NGC 7213). In what follows, and as noted before, we 
collectively refer to this class of galaxy without depleted 
cores as 'Sersic' galaxies, as done by Trujillo et al. (2004) 
who first reported that NGC 4458 and NGC 4478 are 
Sersic galaxies, i.e. they do not display any downward 
departure from the Sersic i? 1 /™ profile at their centres. 
Kormendy et al. (2009) and Hopkins et al. (2009) also 
identify central light excesses over the Sersic function in 
these two galaxies. Except for NGC 4486B (Lauer et 
al. 1996) and NGC 7213, the rms residual scatter is < 

0. 032 mag arcsec -2 for these seven Sersic galaxies. These 
two exceptions are two of only four galaxies, from the 
full sample of 39, with complicated structure (see Ap- 
pendix A Fig. |2"3"|) . As shown in section 6, these seven 
galaxies stand out from the 'corc'-galaxies in a number 
of systematic ways. 

There is some variation in a from galaxy to galaxy 
when using the core-Sersic model. The role of the pa- 
rameter a is to moderate the sharpness of the transition 
between the outer Sersic profile and the inner power- 
law, with higher values corresponding to sharper transi- 
tions and vice versa. The profile of NGC 4291, for exam- 
ple, has a sharper transition and hence requires a larger 
value of a than say NGC 1399. We have set a = 10,5 
and 2 for matching sharp, moderate and broad transi- 
tion regions respectively. The change in the core size, 

1. e. the break radius, while varying a has been closely in- 
spected. The robustness of the break radius is illustrated 
by the plots shown in Fig. [2] which compare the break 
radii obtained using different values of a. The Pearson 
correlation coefficient (r) between the core-Sersic break 
radius Rb,cs( a = 10) and the core-Sersic break radius 
Rb,cs{ a = 5) shown in Fig. Eh> is 1.0. Although the break 
radii from the broad-transition model (i.e. using a = 2) 
has a small amount of scatter in Fig. EK this is due to 
the fact that the majority of our sample galaxies do not 
prefer the broad-transition model. There are, however, a 
small handful of galaxies (NGC 1016, NGC 1399, NGC 
2300, NGC 3379, NGC 4365, NGC 4649, NGC 5419 and 
NGC 5813) whose profiles have a broad transition region 
(see Table EJ. 

4 The Nuker model parameterization (Lauer et al. 2005) reported 
that NGC 4473 and NGC 5576 had break radii of 4.45 and 4.18 
arcseconds, respectively. 




R b cS (a=2) [arcsec] R b cS (a=10) [arcsec] 



Fig. 2. — Robustness of the core-Sersic break radii. Break radii 
from moderate-transition, a = 5, core-Sersic fits plotted against 
a) break radii from broad-transition fits, a = 2 and b) sharp- 
transition fits, a = 10. The Pearson correlation coefficient (r) 
between Rf, c g (a = 5) and Ri, c g (a = 2) is 0.99, while r of Ri, c g 
(a = 5) versus Rb, c S ( a = 10) is 1.0. 

We also estimate the uncertainties of the core-Scrsic 
model parameters by exploring their stability for differ- 
ent a values (i.e. 2, 5 and 10). In agreement with Rich- 
ings et al. (2011) and Trujillo et al. (2004), for some 
galaxies, we notice that 7, n and R e obtained from the 
best a = 5 and 10 fits are slightly different from the ones 
which are obtained using a = 2. 

In addition, we explore the coupling between the core- 
Sersic break radius Rb and the Sersic index n by hold- 
ing the other parameters constant at their best-fit values 
(Tabic 2) and using the x 2 distribution, which we nor- 
malised at the minimum value. In Fig. [3] we show the 
68.3% (ler) confidence limits around the optimal Rb and 
n values for the core galaxies using the A% 2 = 1.0 con- 
tours after marginalizing over the remaining core-Sersic 
parameters. In general, Fig. [3] indicates the absence of 
a coupling between Rb and n, as seen by the small con- 
tours. Overall, we estimate the uncertainties associated 
with 7, Rb, n and R e to be roughly 10%, 10%, 15% and 
20% respectively. We also note that the errors in the 
Sersic parameters could partly be associated with the 
limited radial extent of our data (see section 4.1). 

As noted in passing above, only 4 galaxies are some- 
what poorly represented using our models (NGC 4073, 
NGC 4486B, NGC 6876 and NGC 7213: see Appendix A, 
Fig. EH). NGC 7213 is a Seyfert (Sa) galaxy with a wind- 
ing nuclear dust spiral which can be traced all the way 
to the nucleus (Deo, Crenshaw & Kraemer 2006). It ap- 
pears that this prominent dusty nuclear feature in the 
galaxy profile (see Fig. HJ, which the core-Sersic model 
is not designed to recover, is the likely origin for the 
residual pattern about the model fit to this galaxy, par- 
ticularly in the inner R < 0".3 region. The cD galaxy, 
NGC 4073, has an inner ring (Lauer et al. 2005) over 
the 0".l < R < 0."4 region and we exclude these few 
data points from the fit. NGC 6876 is a dominant el- 
liptical galaxy in the Pavo group with a possible past 
or ongoing merger history (Machacek et al. 2005). The 
residual structures outside the core regions of NGC 4073 
and NGC 6876 that are seen in Appendix A's Fig. [23] 
seem to be associated with the change in ellipticities of 
these galaxies as presented by Lauer et al. (2005). Fi- 
nally, although not apparent from the / and V broad- 
band HST /WFPC2 image, Lauer et al. (1996) remarked 
on the presence of a double optical nucleus from the de- 
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Fig. 3. — Parameter coupling between the core-Sersic break 
radius R\> and Sersic index n for our 'core' galaxy sample. The 
contours show the Ax 2 = 1 boundaries, and their projections onto 
the horizontal and vertical axes give the 68.3% (l<r) confidence in- 
tervals around each galaxies' optimal n and Ri, values, respectively. 
Each contour is generated by holding all the core-Sersic parame- 
ters, other than Ri, and n, fixed at their best-fit values (Table 2). 

convolved WFPC2 image of NGC 4486B. This creates a 
spurious depleted core as noted by Lauer ct al. (1996), 
see also Soria et al. (2006) and Ferrarese et al. (2006). 



4.1. Literature comparison of core-Sersic fits 

We have seven galaxies (NGC 1700, 4291, 4458, 4478, 
5557, 5576 and 5982) in common with Trujillo et al. 
(2004) and 10 galaxies (NGC 4365, 4382, 4406, 4458, 
4472, 4473, 4478, 4486B, 4552 and 4649) in common with 
Ferrarese et al. (2006) and Cote et al. (2006). While Tru- 
jillo et al. (2004) classified NGC 1700 as a Sersic galaxy 
after fitting a profile sampled from R ~ 0".l to 70". 0, 
we tentatively identify a small core within Rb ~ 0".07. 
Apart from NGC 1700, our profile classifications are in 
agreement with Trujillo et al. (2004) and Ferrarese et al. 
(2006). 

The break radii presented in Trujillo et al. (2004) agree 
with our values for all 3 'core' galaxies that we have in 
common. While 4 of the break radii from the 6 'core' 
galaxies in common with Ferrarese et al. (2006) agree 
with our values, NGC 4382 and NGC 4552 are discrepant 
(see Figure [5]). As noted by Ferrarese et al. (2006, their 
section 5.2), the outer disk in the peculiar SO galaxy 
NGC 4382 can, if not modelled as a separate compo- 
nent when the data extends into the domain of the disk, 
bias (high) the Sersic index that would otherwise be as- 
cribed to the bulge of this galaxy. We have therefore 
modelled their extended light profile for NGC 4382 with 
a core-Sersic plus outer exponential disk to show this. 
Our break radius, and Sersic index (see Figure [6]), differ 
from the values presented in Ferrarese et al. (2006) be- 
cause this galaxy's outer disk did indeed bias their analy- 



sis (Figure [7]). Modelling the light profile from Ferrarese 
et al. (2006) for NGC 4552, we recover their fit when 
using their published parameters but find that it can be 
substantially improved upon with a smaller break radii 
and a value of a = 2 (Figures [22] and [7)) . This should be 
compared with Kormendy et al. (2009, their figure 56) 
which reports a core radius in excess of 1 arcsecond from 
their visual inspection. Having accounted for NGC 4382 
and NGC 4552, Fig. 4 reveals an excellent agreement 
between the break radii obtained by us, Trujillo et al. 
(2004) and Ferrarese et al. (2006). Marginal discrepan- 
cies in the break radii may also arise because of the dif- 
ferent filters used in each study. While we use data from 
the F555W filter, Trujillo et al. (2004) primarily used 
data from a F702W filter, and we have taken the F475W 
data, rather than the F850LP data, from Ferrarese et al. 
(2006) as it most closely matches our F555W data. 

When the radial extent of one's data does not ade- 
quately probe the curvature of a spheroid's light pro- 
file, one may not recover the correct Sersic index. Gra- 
ham et al. (2003) noted that truncating a profile from 
^(2-3)i? e to ~ lR e changes the fit parameters by up to 
5%. Concerned about this, as our data only extends to 
10", we have additionally compared our Sersic indices 
with those from Trujillo et al. (2004) and Ferrarese et al. 
(2006) which had a larger radial extent. First, we note 
that the Sersic index for NGC 4458 from Trujillo et al. 
(2004) was biased (high) by the presence of a nuclear disk 
that was not separately modelled as we have done here. 
They noted that they were not confident in their analy- 
sis of this galaxy and as such we have excluded this one 
galaxy from Trujillo et al. (2004) in our Figure |5J With 
the exception of NGC 4552 and NGC 4382, mentioned 
in the preceding paragraph, the agreement between the 
Sersic indices is good to within 50 per cent or better, 
which in fair agreement with the la uncertainty range of 
±36 percent found by Allen et al. (2006, their figure 15). 
However, it should also be noted that due to ellipticity 
gradients, the major-, minor- and geometric-mean axis 
do not have the same Sersic index (Ferrari et al. 2004). 
Their values can disagree by up to a factor of ^2 (e.g. 
Caon et al. 1993, their figure 4). Consequently, some of 
the scatter seen in Fig.[6l which compares our major-axis 
Sersic indices with the geometric-mean axis values from 
Ferrarese et al. (2006) is because of this. 

We have also been able to include a comparison of 
break radii and Sersic indices, in Figs [5] and [6] respec- 
tively, for 8 core galaxies that we have in common with 
Richings ct al. (2011). In general, the agreement is good, 
although there are three somewhat discrepant points. 
For NGC 5982, we suspect that Richings et al. may have 
missed the core with their Sersic n = 2 fit to this large 
elliptical galaxy with a = 239 km s _1 and My = —22 
mag. We are also inclined to prefer our n = 4 core-Sersic 
fit to NGC 5813 rather than the n = 9 core-Sersic fit of 
Richings et al. which results in a larger apparent core. 
However we suspect that the dust ring in NGC 3607, 
although subtracted, may have still interfered with the 
optical light-profile from Lauer et al. (2005). As such, 
we feel that our break radius may be over-estimated for 
this one galaxy, explaining its outlying nature in the cen- 
tral mass deficit versus black hole mass diagram (see sec- 
tion [73]). 

Recently, Dhar & Williams (2011) noted that surface 
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Fig. 4. — Optical HST images of galaxy centres. Each galaxy was observed using the high-resolution (0".046/pixel) PC1/WFPC2 
camera. 
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Fig. 5. — Comparison of the core-Sersic model's break radii from 
this study (ijj,: Table 2, V-band) and previous break radii, Rf, }P rev) 
from i) Trujillo et al. (2004, if-band, filled stars), ii) Ferrarese 
ct al. (2006, g-band, open stars) and iii) Richings et al. (2011, 
various bands, filled circles). We have converted the geometric 
mean radii from Ferrarese et al. (2006) into semi-major axis radii 
using their cllipticity values. A Representative error bar is shown 
at the bottom of the panel. 

brightness profiles of elliptical galaxies can be modelled 
well using the multi-component DW-function (Dhar <fc 
Williams 20f0), with smaller rms residuals than the 
Sersic, Core-Sersic and Nuker models. Although we fit 
the profiles (from Lauer et al. 2005) with limited radial 
ranges, for galaxies that we have in common with Dhar 
& Williams (2011) — NGC 4365, NGC 4382, NGC 4406, 
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Fig. 6. — Comparison of our major-axis Sersic indices, n, with 
previously published values (n prev ) derived using a greater radial 
extent. The filled stars and circles are major-axis Sersic indices 
from Trujillo et al. (2004) and Richings et al. (2011) respectively 
while the open stars are geometric-mean axis Sersic indices from 
Ferrarese et al. (2006). 

NGC 4458, NGC 4472, NGC 4473, NGC 4478, NGC 4552 
and NGC 4649 — our rms residuals arc smaller than 
those obtained with the DW-function (Dhar & Williams 
2011, their Table 1) except for NGC 4458 (Appendix 
A Fig. 

5. BREAK RADII MEASUREMENTS 

Given that the core-Sersic model and the Nuker model 
yield different core sizes, we have included an extended 
section on the measurement of the break radii via differ- 
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TABLE 2 

Structural parameters 



Galaxy 


^b,V 


Rb 


Rb 


7 


a 


n 


Profile Type 


m pt ,v 


fio v h Notes 




_9 

mag arcsec 


(arcsec) 


(PC) 










(mag) 


mag arcsec (arcsec) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) (11) (12) 


NGC 0507 


16.45 


0.42 


130 


0.09 


5 


3.2 


c-S 






NGC 0584 


14.57 


0.19 


18 


0.46 


5 


6.4 


c-S 






NGC 0741 


17.52 


0.96 


338 


0.19 


5 


7.4 


c-S 


22.4 




NGC 1016 


17.21 


0.68 


289 


0.20 


2 


5.2 


c-S 






NGC 1374+ 


— 


— 


— 


— 


— 


2.8 


s 




14.3 0.17 


NGC 1399+ 


16.29 


2.09 


196 


0.11 


2 


4.0 


c-S 






NGC 1700 


13.64 


0.07 


18 


0.27 


5 


8.9 


c-S 






NGC 2300 


16.85 


0.98 


122 


0.11 


2 


1.7 


c-S 






NGC 3379 


15.66 


1.03 


51 


0.19 


2 


4.0 


c-S 






NGC 3607 


16.42 


1.52 


164 


0.29 


5 


5.0 


c-S 




nuclear dust lanes 


NGC 3608 


15.10 


0.21 


23 


0.28 


5 


5.1 


c-S 






NGC 3640 


14.72 


0.04 


5 


-0.02 


5 


3.9 


c-S 






NGC 3706 


14.17 


0.11 


21 


-0.02 


10 


6.3 


c-S 




ring of stars (0".06-0".4) 


NGC 3842 


17.42 


0.72 


320 


0.19 


5 


6.1 


c-S 






NGC 4073 


16.47 


0.24 


99 


-0.05 


10 


8.1 


c-S 




ring of stars (0".l-0".4) 


NGC 4278 


15.79 


0.75 


56 


0.20 


5 


1.1 


c-S 


19.4 


NGC 4291 


15.22 


0.35 


11 


0.10 


5 


5.3 


c-S 






NGC 4365 


16.56 


1.40 


135 


0.04 


2 


1.1 


c-S 


20.1 




NGC 4382 


15.04 


0.32 


28 


0.08 


5 


3.7 


c-S 






NGC 4406 


16.01 


0.87 


70 


0.01 


5 


8,1 


c-S 






NGC 4458 


— 


— 


— 


— 


— 


3.1 


S 


17.0 


14.4 0.25 


NGC 4472 


16.44 


1.68 


129 


0.00 


2 


1.3 


c-S 


22.2 




NGC 4473 


— 


— 


— 


— 


— 


2.1 


s 






NGC 4478 


— 


— 


— 


— 


— 


2.7 


s 


20.1 


15.6 0.41 


NGC 4486B 


— 


— 


— 


— 


— 


3.0 


s 




double optical nuclei 


NGC 4552 


15.00 


0.35 


25 


0.01 


10 


3.9 


c-S 


20.5 




NGC 4589 


15.34 


0.21 


22 


0.30 


5 


5.5 


c-S 






NGC 4649 


16.92 


3.23 


256 


0.18 


2 


3.6 


c-S 






NGC 5061 


14.06 


0.21 


33 


0.13 


5 


7.0 


c-S 






NGC 5419 


17.60 


1.67 


485 


-0.05 


2 


6.6 


c-S 


19.9 




NGC 5557 


15.31 


0.16 


35 


0.17 


5 


3.1 


c-S 






NGC 5576 












3.5 


S 






NGC 5813 


16.15 


0.42 


64 


-0.09 


2 


1.1 


c-S 






NGC 5982 


15.45 


0.21 


12 


0.08 


5 


3.3 


c-S 






NGC 6849 


16.72 


0.22 


84 


0.20 


5 


6.3 


c-S 






NGC 6876 


16.98 


0.45 


119 


-0.01 


10 


6,1 


c-S 




double optical nuclei 


NGC 7213 












1.5 


S 


16.6 


12.6 0.04 


NGC 7619 


15.78 


0.35 


87 


0.12 


5 


1.1 


c-S 






NGC 7785 


14.94 


0.05 


11 


0.16 


10 


5.3 


c-S 







Notes. — Structural parameters from fits to the V-band major-axis surface brightness profiles (Appendix A Figs. [221 and I23jl . The superscript + 
is used to indicate that an F606W surface brightness profile is used, rather than an F555W surface brightness profile. Col. (1) Galaxy name. 
Col. (2)-(7) Best-fit parameters from the core-Sersic model, equation (2). Col. (8) Indicates the profile classification where c-S = core galaxy 
described by the core-Sersic model, and S= Sersic galaxy described by the Sersic model. Col. (9) Point source apparent magnitude. Col. (10) 
Nuclear disk central surface brightness. Col. (11) Nuclear disk scale length. Col. (12) Description of inner additional light components. A "?" 
indicates a tentative classification. 
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5.1. The Nuker model 



As noted previously, the Nuker team (Lauer et al. 1995, 
2005) used a 5-parameter double power-law model for fit- 
ting the inner radial surface brightness profiles of galax- 
ies. Dubbed the "Nuker law" , it can be written as 



I(R) = 7 6 2^-tV q 



R 
Rb 



R 
~R, 



(7-/3)/" 



(4) 

where is the intensity at the break radius i?t, . The neg- 
ative logarithmic slopes for the inner and outer power- 
law regions are denoted by 7 and /3, respectively, while 
a controls the sharpness of the transition. Rf, represents 
both the radius of maximum curvature of this model, and 
the location where the local gradient of the model equals 
— (/3 + 7)/2. One can readily appreciate how fitting this 
model to a larger radial range, and thus to an increas- 
ingly steeper outermost region of what are curved sur- 
face brightness profiles, results in larger values of /3 and 



thus larger values of Rb- When this occurs, the value of 
a is also reduced, to accommodate an (artificially) in- 
creasingly broad transition. Having excluded additional 
nuclear components, Lauer et al. (2005) tabulated their 
best fitting Nuker model parameters for the galaxies used 
in this study, see also Lauer et al. (2007a, b), and we refer 
to those values in some of the following figures. 

In their study of the nuclear regions of early type galax- 
ies, Rest et al. (2001) noted that for small values of a 
(i.e. broad transitions), the Nuker model's parameter 7 is 
rather a representation of the slope of the brightness pro- 
file at radii much smaller than the image resolution limit. 
They therefore introduced another parameter, 7', that 
was the negative logarithmic slope of the Nuker model 
at R = 0".l, and which was adopted by Lauer et al. 
(2005, 2007b). Graham et al. (2003) had however noted 
that because this local logarithmic slope 7' (at R = 0".l) 
is a distance dependent quantity, galaxies with identical 
surface brightness profiles observed at different distances 
will have different 7' values. That is, this quantity is not 
a physically robust or meaningful quantity to use. 
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Fig. 7. — Left panels: Core-Sersic models presented in Fer- 
rarese et al. (2006). Right panels: Our modelling of their pub- 
lished profiles, including an outer exponential disk (dotted line) 
for NGC 4382. Note: while the innermost data points of the flat 
cores are slightly affected by the PSF, which we have not accounted 
for here and thus things look worse for the Ferrarese ct al. fit than 
they are at small radii, this is not the reason for the different Sersic 
indices and break radii. For the four other core-Sersic galaxies that 
we have in common with Ferrarese et al. (2006), we agree on the 
location of the break radii (Figure [5}. 

In Fig. |8] we show how the 7 values from the core-Sersic 
model fit to 32 galaxies (having depleted cores relative to 
the outer Sersic profile) compare with the 7 (Pearson's 
r = 0.53) and 7' (Pearson's r = 0.87) values obtained 
from the Nuker model. While the latter correlation co- 
efficient is high, it is pointed out that this does not indi- 
cate a slope of unity, which would be required if the two 
parameters were equivalent. Instead, for positive values 
of 7 c s, the plots reveal that the (Nuker model)-derived 
values of 7^ u k an d 7Nuk are generally smaller than the 
core-Sersic 7 values. 

5.2. Model independent approaches 
5.2.1. — d 2 log //dlog R 2 = maximum 

In an effort to obtain a nonparamctric estimate of the 
break radii, we attempt to locate the radius correspond- 
ing to the maximum of the second derivative of the ob- 
served intensity profile in logarithmic coordinates, inde- 
pendent of any model or any smoothing or alteration 
of the data. In principle this model-independent ra- 
dius should mark the break radius. In practice, how- 
ever, Fig. [9] reveals that there exists several compara- 
ble maxima over an extended radial range due to the 
sensitivity of this approach to the noise in the profile 
data. Therefore, we were unable to use this technique 
to acquire accurate break radii. We do however show in 
Fig. ([TU1 Right) that the core-Sersic model's break radius 
corresponds to the radius where the second logarithmic 
derivative of this model has its maximum value. 

In passing we note that Lauer et al. (2007a, their fig- 
ure 17) and Lauer et al. (2007b, their figure 16) had re- 
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Fig. 8. — Comparison of the core-Sersic model's inner negative 
logarithmic slopes (7cSi Table 2) with a) the Nuker model's inner 
negative logarithmic slopes (7jv u fc, from Lauer et al. 2005), and 
b) the Nuker model's local negative logarithmic profile slope at 
the instrument resolution limit {~l' Nuk , from Lauer et al. 2005). 
Representative error bars are shown at the bottom of each panel. 

ported a strong consistency between their Nuker model 
break radii Rb,Nuk and the location of the maximum of 
the second derivative of the intensity profile, in logarith- 
mic coordinates, acquired from a model-independent ap- 
proach. This claim was however surprising because it is 
well established and understood why Rb,Nuk varies con- 
siderably as the fitted radial extent of a galaxy's surface 
brightness profile is varied (e.g. Graham et al. 2003, their 
figures 2-4). As such, while the varying Rb,Nuk always 
corresponds to the radius where the fitted Nuker model 
has the maximum of its second derivative, this typically 
will not correspond to the radius where the actual data 
has the maximum in its second derivative. Using the 
same data as Lauer et al. (2007a, b), we cannot repro- 
duce their result. In Fig. 8, we have shown that the sec- 
ond derivative of the observed galaxy brightness profiles 
has a large point-to-point variation. It is not clear how 
the result of Lauer et al. (2007a, b) could be obtained 
without smoothing the data. Furthermore, we continue 
to reiterate the result of Graham ct al. (2003), that the 
Nuker model fits, especially the break radius, are depen- 
dent on the chosen fitting range. In spite of this, that 
work was the sole basis for their rejection of the concerns 
about the Nuker model raised by Graham et al. (2003). 

5.2.2. -dlog 7/ dlog R = 7' = 1/2 

In a continued effort to better measure the sizes of 
partially depleted galaxy cores, this section investigates 
an alternative model-independent radius that has been 
used in the literature. This investigation is important 
if we are to accurately quantify the extent of damage 
caused by coalescing SMBHs at the centres of galaxies. 

Most 'core' galaxies have negative, logarithmic, in- 
ner profile slopes less steep than 0.3-0.5 (Glass ct al. 
2011, and references therein), before they transition to 
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brightness profile (in theV-band) plotted against semi-major axis 
(see section [5.2. It . 

an outer Sersic profile with a slope typically steeper than 
0.5 (Lauer et al. 2005) 5 . The negative logarithmic slope 
beyond the core, in the Sersic portion of the profile, varies 
as {b n /n)(R/R c ) 1 / n . Graham et al. (2003, their Figure 6) 
have revealed that this slope is steeper than 0.5 for Sersic 
models with n > 3-4 once beyond 1% of the effective ra- 
dius R e . Given that Cote et al. (2007) have reported 

that cores extend to ~0.02to;oi 5 -ffe, we can appreciate 
why the slopes on the outer side of cores are steeper 
than 0.5. As such, the radius where the negative loga- 
rithmic slope of the underlying galaxy surface brightness 
profile 7' equals 1/2 (Carollo et al. 1997) can be used 
to find the transition radius between the inner core and 
the outer Sersic profile, and thus approximate the break 
radius of the core-Sersic model. 

In finding the radius where 7' = 1/2, when present 
we avoided data points affected by central light excesses 
before applying this technique (Fig. [TU1 Left). We visu- 
ally inspected individual profiles to double-check these 
objective break radii estimates and for (only) five galax- 

As discussed by Graham &; Guzman (2003), pure 'Sersic' 
galaxies can also have inner profile slopes less steep than 0.5, and 
whether or not one measures this (underneath any additional nu- 
clear components) simply depends on how small one's R/R c spatial 
resolution is (see Graham et al. 2003, their Figure 6). 
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Fig. 10. — Left: Typical examples of the negative logarith- 
mic slopes of the surface brightness profile (in the l/-band) plotted 
against semi-major axis (see section [5.2.21 1. Right: Typical exam- 
ples of the second logarithmic derivative of the fitted core-Sersic 
model (in the V-band) plotted against semi-major axis. Arrows 
indicate the break radii of the galaxies from the fitted core-Sersic 
model (Table [2}. 

ies (NGC 3607, NGC 3640, NGC 4406, NGC 4589 and 
NGC 5557) the data were too noisy; as such we excluded 
them from the comparison of break radii in the follow- 
ing subsection. The radii where 7' = 1/2 are plotted in 
Fig. [UJ and discussed in the following subsection. 

Before proceeding, we note that Lauer et al. (2007a) 
had also remarked that Carollo et al. (1997) had advo- 
cated use of Rryi, with 7' = 1/2, as a core scale param- 



eter. They wrote that "Since R 



7 '=l/2 



is generally well 



interior to Rb,Nuk, it is not meant to describe the ac- 
tual complete extent of the core; it is just a convenient 
representative scale." To reduce potential misinterpre- 
tation of this comment, it is important to note that the 
Nukcr model's Rb,Nuk was also never meant to describe 
the complete extent of the core. While the published 
Nukcr model break radii occur at large radii where the 
slope is steeper than —0.5 (see Fig. [Tib), the actual outer 
edge of the Nuker model's transition region occurs con- 
siderably further out than this, greater than both the 
Nuker model break radii and the outer edge of the core- 
Sersic model. Moreover, due to the curved nature of the 
outer Sersic profile, the outer edge of the Nuker model 
increases as the fitted radial extent is increased. 

5.3. The core-Sersic model 

Fig. [TTb reveals a strong correlation (Pearson's r = 
0.96) between the break radii estimated from the model- 
dependent (core-Sersic) and the mo del- independent 
(R-fi =1/2) quantitative analysis. One can similarly show 
that the negative logarithmic slope of the fitted corc- 
Sersic models has a value of 0.5 at a radius very close 
to the core-Sersic model's break radius. This important 
result has not been noted before. Turning things around, 
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Fig. 11. — Comparison of R i = 1/2 and a) the break radius 
from the core-Sersic model (R b c g , see Table 2) and b) the Nuker 
model (Rf, jvufci Lauer et al. 2005). When present, central light 
excesses are avoided when determining R 1 r— 1 / 2 - NGC 3607, NGC 
3640, NGC 4406, NGC 4589 and NGC 5557 are excluded from this 
analysis due to considerable noise in their inner data. A Represen- 
tative error bar is shown at the bottom of panel a. 
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Fig. 12. — Comparison of the published Nuker break radii (Lauer 
et al. 2005) and the core-Sersic break radii (see Table 2). A repre- 
sentative error bar is shown at the bottom of the panel. 
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Fig. 13. — The ratio of the Nuker break radius (Rb,Nuk> Lauer 
et al. 2005) to the core-Sersic break radius (iJt, jC s, see Table 2) for 
core galaxies as a function of a) Rb }C S an( i b) the core-Sersic model 
negative inner logarithmic slope f a s- 



this agreement endorses the use of R-y'=i/2 when the data 
are not too noisy, albeit with the caveat that even galax- 
ies without partially depleted cores may still have a re- 
solvable radius where 7' = 1/2. The core-Sersic model 
is therefore still recommended for identifying and quan- 
tifying cores. 

For comparison, Figs. ITTb and LT21 show that the mea- 
surement of the Nuker model's break radii presented in 
Lauer et al. (2007a, b) against the model-independent 
break radii -ff 7 < = i/2 (r = 0.91) and the core-Sersic break 
radii (r = 0.90). Graham et al. (2003) explained that 
increasing the fitted radial range of the Nuker model, 
into the curved profile beyond the core, will increase the 
slope of the Nuker model's outer power-law (3 and thus 
result in the Nuker model's break radius correspond- 
ing to the location where the model's slope is the aver- 
age of the inner and outer power-law slopes 7 and /3 — 
marching out to larger radii. For this reason the Nuker 
model's break radius can be pulled out beyond the actual 
transition radius between the core and the outer Sersic 
profile. On average, the Nuker break radii are ~2 times 
bigger than the core-Sersic break radii (Fig [13]). While 
Lauer et al. (2007a, their Appendix C) refuted this crit- 
icism of the Nuker model, they simultaneously reiter- 
ated this very problem and subsequently used i? 7 ' = i/2- 
However they did not show how i?y =1 / 2 compared with 
the Nuker model break radii, which is done here for the 
present galaxy sample that have 'real' cores (Fig ITTb ) . 
The above problem with the Nuker model was, in part, 
the motivation for the core-Sersic model. Fig. [TJ2 re- 
veals that there is no correlation between the ratio of the 
Nuker break radius (Rb,Nuk) and the core-Sersic break 
radius (Rb, c s) with either a) Rb, c s (Pearson coefficient 
r = —0.29, Fig [T3a ) or b) the core-Sersic negative in- 
ner logarithmic slope 7 c s (Pearson coefficient r= -0.12, 

Fig. nab). 

From Table [JJ we see that the break radii are smaller 
than 0.5 kpc. It would appear that the 500 pc resolu- 
tion models by Martizzi et al. (2012, their Figure 7) may 
have 'over-cooked' core-formation in their simulations of 
galaxies with cores up to ~8-10 kpc in size. Similarly, 
the large ~3 kpc cores created by Goerdt et al. (2010) 
are not observed in real galaxies. 

6. IDENTIFICATION OF DEPLETED CORES, AND THEIR 

SLOPES 

Initially, Kormcndy et al. (1994) and Lauer et al. 
(1995) identified cores if the inner slope of the Nuker 
model was less than 0.3. Kormendy (1999) subsequently 
relaxed this criteria to read "galaxies that show a break 
from steep outer profiles to shallow inner profiles" , with 
the "outer profile" modelled by the outer power-law of 
the Nuker model. However this definition of a core, and 
modelling of the stellar distribution, resulted in a dis- 
connection with the curved outer Sersic profile that was 
known to exist (e.g. Caon et al. 1993). Graham et al. 
(2003) therefore advocated that cores be identified and 
defined as a central stellar deficit relative to the outer 
Sersic profile. Kormendy et al. (2009) quoted and par- 
tially embraced this new definition but opted to identify 
by eye the region to fit the Sersic model and thus the 
onset of the break radius, rather than using the core- 
Sersic model in an objective analysis. Their visual core 
classification agrees with our core identification, however, 
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their approach resulted in break radii notably larger than 
those given by the core-Sersic model (cf. Graham 2004; 
Ferrarese et al. 2006) and thus, given the results in the 
previous section, their break radii do not agree with the 
model-independent measurement of where the light pro- 
file has reached a steep outer profile with negative loga- 
rithmic slope equal to 0.5. This arose in part due to the 
lack of an infinitely sharp transition region. These larger 
break radii also result in larger estimates of the central 
mass deficit. Using a model- independent technique based 
on average profiles, Hopkins & Herquist (2010) reported 
an upper limit to the mass deficit of 2-4 times the central 
black hole mass, in agreement with the analysis presented 
in Graham (2004) using the core-Sersic model. This is 
in contrast to the values of 10-20 times the central black 
hole mass reported by Kormendy & Bender (2009). 

Although all of the 39 galaxies in our sample were tab- 
ulated as 'core' galaxies by the published analysis using 
the Nuker model (Lauer et al. 2005, their Table 4) 6 , in 
section 2] we effectively reclassified seven as Sersic galax- 
ies without partially depleted cores. This concern over 
misidentification in lower luminosity spheroids was first 
highlighted by Graham et al. (2003), and the discrepancy 
is fundamentally due to the inclusion of Serisc galaxies 
with low Sersic index n (and thus a shallow inner profile 
slope) that have no depleted core relative to the Sersic 
profile which describes the outer galaxy light distribu- 
tion. The typical value of the Sersic shape index (n) for 
the seven Sersic galaxies is w 3 (Fig. [T4|) . 

In Figure [15] we have plotted the central surface bright- 
ness, ^o,Vj °f the galaxies or bulges obtained from the 
fitted models against their T^-band absolute magnitude 
My. As with the My-n diagram (Fig. [14]) . it is immedi- 
ately apparent that the seven non-core galaxies are not 
some random sample from the 39 galaxies. They reside in 
a region of the My-fx^y diagram known to be occupied 
by galaxies without partially depleted cores. 

Prior to Lauer et al. (2005), Trujillo et al. (2004) 
had already revealed that two of these seven galaxies 
(NGC 4458 and NGC 4478) had no partially-depleted 
core. Ferrarese et al. (2006) additionally identified that 
NGC 4473 did not have a depleted core, and Kormendy 
et al. (2009) subsequently acknowledged that none of 
these three galaxies have depleted cores but instead have 
an excess of nuclear flux. In addition to the double- 
nucleus galaxy NGC 4486B which is known to have a 
false core (Lauer et al. 1996), we find that NGC 1374, 
NGC 5576 and NGC 7213 also do not have cores depleted 
of stars relative to their host spheroid's (outer) Sersic 
profile. NGC 7213 is a rather faint spiral galaxy with 
Mb ~ —19.5 mag and is thus not expected to have a de- 
pleted core like luminous, boxy spheroids do. NGC 1374 
is also not a luminous galaxy; it too has Mb ~ —19.5 
mag, and NGC 5576 is only 0.7 mag brighter. 

Besides their distribution in Figs. [14] and [15] it is of 
interest to examine if there are additional characteristics 
among these seven galaxies which have cores according to 
the Nuker model (Lauer et al. 2005) but do not have par- 
tially depleted cores relative to their outer Sersic profile. 
From the full sample of 39 galaxies, nine have velocity 

6 Although tabulated as having a 'core' profile in Lauer et al. 
(2005), those authors are aware that NGC 4486B is not a 'real' 
core galaxy (Lauer et al. 1996). 
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Fig. 14. — Absolute l/-band galaxy magnitude from Table 1 
(converted to bulge magnitude for the lenticular and spiral galax- 
ies) plotted as a function of the Sersic index, n, that quantifies 
the shape of the underlying galaxy or bulge, major-axis light pro- 
file. The line Mg = — 9.41og(n) — 14.3 is taken from Graham & 
Guzman (2003, their figure 10) and adjusted here to the V-band 
using B — V = 0.9 (Fukugita, Shimasaku & Ichikawa 1995). Filled 
circles are the core galaxies; filled triangles are the Sersic galaxies; 
bulges are circled. The Pearson coefficient between M.y and n is 
r= —0.34. A representative error bar is shown at the bottom of 
the panel, but the lcr uncertainty on the bulge magnitudes for disk 
galaxies is ~ 0.75 mag. 

dispersions a < 183 km s~ x according to HyperLeda, and 
seven of these nine galaxies are the above Sersic galaxies. 
Of the remaining two galaxies, NGC 3640 has a = 182 
km s _1 and a very small, questionable depleted core, and 
while NGC 4382 has a = 179 km s _1 according to Hyper- 
Leda's mean value however the most recent measurement 
of a = 205 ±8 km s _1 (Bcrnardi et al. 2002) is supportive 
of a 'real' core. 

Aside from the above mentioned confirmation by other 
authors, there are good reasons to suspect that the Nuker 
model struggles to identify cores which have been de- 
pleted relative to the inward extrapolation of the outer 
Sersic profile in galaxies and bulges fainter than My w 
-21 mag (see Figure Il4]) . While Lauer et al. (2007a) 
used their Nuker model parameters to determine where 
the fitted Nuker model has a slope of —1/2 which we 
consider to be a better measurement of the transition 
from a shallow core to a steep outer profile — one is left 
with the problem of not knowing which galaxies actually 
have depleted cores, that consistent with non-parametric 
identification methods, relative to the smooth inward ex- 
trapolation of the outer Sersic profile. Figure [141 reveals 
that knowledge of this Sersic index can help with this 
diagnosis, but it obviously requires remodelling the light 
profiles. An additional criteria, which could be better 
quantified with a larger galaxy sample, is that alleged 
cores in galaxies with a < 183 km s" 1 are likely not to 
correspond to a real deficit of stars. 

Although our sample is not suitable for testing the bi- 
modality, or lack thereof (Glass et al. 2011, and references 
therein) , in the distribution of host galaxy central surface 
brightness slopes against absolute magnitude (which, for 
the reason we have just seen, can not be used as a diag- 
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Fig. 15. — Absolute V-band magnitude of the galaxy (or the 
bulge of the disk galaxies) plotted against the core-Sersic model's 
inner V-band surface brightness fio at the instrument resolution 
limit R = 0".0455 (excluding NGC 1374 and NGC 1399 for which 
there were no F555W profiles). Filled stars and triangles are 
core-Sersic and Sersic galaxies respectively. Galaxies with addi- 
tional nuclear light have the host galaxy plus central excess flux, 
at R = 0".0455, shown by the open symbols. Solid lines connect 
the central brightness values of the nucleated galaxies and their 
underlying host galaxies. The dotted line is taken from Graham Sz 
Guzman (2003, their figure 9) and adjusted here to the V-band us- 
ing B — V = 0.9. The dashed line, My = -1.09(/i - 16.0) - 22.60, 
shows our least-squares fit to the My and /io relation for galax- 
ies with 'real' depleted cores, and having a Pearson's coefficient 
r = —0.44. A representative error bar is shown at the bottom of 
the panel, but the la uncertainty on the bulge magnitudes for disk 
galaxies is ~ 0.75 mag. 

nostic of core identification), we do find that 7 < 0.3 for 
all (core-Sersic)-identificd core galaxies, with the possible 
exception of only NGC 584 (where 7 = 0.46, see Table 2, 
and Fig. I16j). The low-luminosity (low-n) Sersic galaxies 
without cores can have inner slopes ranging from to 0.5 
and steeper depending on the radius where one samples 
the underlying host spheroid's Sersic profile ( Graham & 
Driver 2005, their section 2.4). Fig. \W\ plots the relation 
between 7 and the V-band absolute magnitude My for 
the 'true' core galaxies given in Table 2. For the first 
time using the core-Sersic model, we find that a number 
of bright (high n) core galaxies exhibit values of 7 less 
than (see also Lauer et al. 2005, and references therein 
for galaxies with Nuker model values of 7 < 0). Earlier 
studies with the core-Sersic model had been confined to 
reporting 7 > 0. Kandrup ct al. (2003) have discussed 
how black hole binaries may couple with the innermost 
stars and transport them to a larger radius, resulting in 
such surface brightness profiles. 

7. STRUCTURAL PARAMETER RELATIONS 

Having selected a set of cores based on the core-Sersic 
model, we intend to obtain radial profiles over a greater 
radial extent for those 32 galaxies in a future paper. 
Fig. [TJ] illustrates the linear correlation between the ab- 
solute galaxy magnitude (bulge magnitude in the case of 
disk galaxies) and the light profile shape, n. As noted 
in section 2, the bulge magnitude for the one spiral (Sa) 
and six lenticular (SO) galaxies in our sample are ob- 
tained using the representative bulge-to-disk flux ratio 
given in Graham & Worley (2008, their Table 5). Al- 
though our galaxies arc limited in number and range of 
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Fig. 16. — Comparison of the core-Sersic model's negative inner 
logarithmic slope 7 with the V-band absolute galaxy magnitude 
(bulge magnitude for disc galaxies). Bulges in disk galaxies are 
circled. A representative error bar is shown at the top of the panel, 
but the lc uncertainty on the bulge magnitudes for disk galaxies 
is ~ 0.75 mag. 

absolute magnitude (-18.77 > My > -23.45 mag), the 
overall distribution in the L—n diagram is in agreement 
with that from Caon et al. (1993), Graham & Guzman 
(2003; their Fig. 10) and Ferrarese et al. (2006). Given 
that the Sersic index n is derived from 10" profiles, we do 
not consider it to be as accurate as possible. Nonethe- 
less, on average, the galaxy ensemble adhere to the es- 
tablished L-n relation, and are in fair agreement with the 
indices derived from fits to a larger radial extent (see sec- 
tion 4.1). This gives us some confidence that our Sersic 
parameters (which are only used once in this paper to 
derive a 'ballpark' result in section \7.3\i are not too far 
off from the correct values. 

7.1. Rb~L, Rb-cr, Rb-^b, f-b-L and [ib-o relations 

Building on earlier works which explored the connec- 
tion between early-type galaxy dynamics and isophotal 
shape (e.g. Davies et al. 1983; Nieto & Bender 1989; Ni- 
eto et al. 1991), Faber et al. (1997) highlighted associa- 
tions between the core structure and global galaxy prop- 
erties. Bright ellipticals with boxy isophotes, slow rota- 
tion and pressure supported dynamics arc 'core' galax- 
ies, while fainter elliptical galaxies with elliptical or 
'discy' isophotes and often rotational support are 'Sersic' 
('power-law') galaxies. They went further to argue that 
the presence of a 'core' is a better predictor of the slow 
rotation or boxiness than the galaxy absolute magnitude, 
although we have just learned that some of the Nuker- 
derived cores are not consistent with core-Sersic cores or 
visual classification, i.e. they have been claimed to exist 
in galaxies with no depleted core. It would be remiss if we 
did not use our refined core-Sersic parameters to derive 
updated scaling relations for galaxies that have central 
stellar deficits relative to their outer Sersic profile. That 
is, we exclude those galaxies having no depleted cores. 

Figs.flTk andfTTb display the relation between the core- 
Sersic break radius (Table 2) and the published Nuker 
'cusp radius' i? 7 > =1 / 2 (Lauer et al. 2007a), collectively 
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Fig. 17. — The core-Sersic break radius and the published Nuker 'cusp radius' ily—j/a (Lauer et al. 2007a), collectively Ri, (in pc), are 
plotted as a function of a) absolute V-band spheroid magnitude (Table [TJ, b) central velocity dispersion a (Table [TJ and c) the associated 
V-haad surface brightness, collectively referred to as fj, b here (excluding NGC 1399 for which there was no F555W profile). Filled stars are 
the core galaxies from this study fitted with the core-Sersic model, while open circles are these same galaxies with and fn, (7' = 1/2) data 
from Lauer et al. (2007a). Triangles show the location of the seven alleged 'core' galaxies which are reclassified here as Sersic galaxies and 
thus have no core-Sersic break radii. The solid lines are least-squares fits to the corc-Serisc data, while the dashed lines are the published 
Nuker relations after they excluded galaxies with My > —21 mag (Lauer et al. 2007a, their Eqs. 13, 14 and 17). Pearson correlation 
coefficients, r, (and representative error bars), for our data, are shown at the top (bottom) of each panel. For the disk galaxies, the la 
uncertainty on the bulge magnitudes is ~ 0.75 mag. 



Rb, and a) the U-band absolute magnitude My listed in 
Table 1, and b) the central velocity dispersion a (Table 
1). Using the ordinary least-squares (OLS) bisector re- 
gression (Feigelson & Babu 1992), a fit to the core-Sersic 
Rb and a yields 
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and further application of the bisector regression gives 
the relation between the core-Sersic Rb and My as 



O 

CD 

m 
o 

M 

crJ 



15 



log 



18 



(-0.58 ± 0.09) {M v + 22)+ (1.90 ±0.10). 



(6) 

Similar trends between break radius and galaxy magni- 
tude were also seen in Faber et al. (1997, their figure 4), 
Ravindranath et al. (2001, their figure 5a, b), Laine et 
al. (2003, their figure 9), Trujillo et al. (2004, their fig- 
ure 9), de Ruiter et al. (2005, their figure 8) and Lauer 
et al. (2007a, their figure 19, bottom panel). The three 
outliers in our R b -M v relation are: NGC 2300 and NGC 
3607 — lenticular (SO) galaxies with a big core for their 
bulge brightnesses — and NGC 3640, a galaxy known for 
its morphological peculiarity which probably signals an 
on going or a recent merger (Michard & Prugnicl 2004). 
Shown in Fig. 117b is Rb (core-Sersic break radius from 
Table 2 plus Nuker 'cusp radius' from Lauer et al. 2007a) 
as a function of fib (the F-band surface brightness at 
the core-Sersic break radius and the Nuker 'cusp radius', 
from Table 2 and Lauer et al. 2007a respectively). The 
bisector fit to the core-Sersic Rb and fib gives 
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Fig. 18. — A compilation of the core-Sersic fits to the major-axis 
surface brightness profiles of all core galaxies from Table 2 except 
for NGC 1399 for which there was only an F606W profile . St ars 
indicate the break radius of the individual profiles (cf. Fig. I17fc ). 

Equations 5, 6 and 8 update the Rb-cr, Rb-L and fib- 
Rb relations presented in Lauer et al. (2007a, their equa- 
tions 13, 14 and 17 respectively). Fig. [TBI and Fig. [i~7b 
confirm, but re-define the tight correlation between the 
core-Sersic core brightness fib and the core-Sersic core ra- 
dius R b (Faber et al. 7 1997; de Ruiter et al. 2005, their 
figure 7; Lauer et al. 2007a, their figure 6). Apparent in 
Fig. [18] is a roughly universal profile beyond the core of 
'core' galaxies, out to ~1 kpc, which explains the tight 

7 Faber et al. (1997) ascribed the tight correlation among the 
central properties of early-type galaxies to the presence of a "core 
fundamental plane" (log rj, fif,, log a), which is analogous to the 
global Fundamental Plane — log r e , /i e , log a (e.g. Djorgovski & 
Davis 1987). 
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Fig. 19. — The core-Sersic model's break surface brightness and the published Nuker's 'cusp' surface brightness, collectively fn, ( V-band), 
are plotted against a) the absolute l/-band magnitude of the galaxy (or the bulge for disc galaxies) and b) the central velocity dispersion 
a (see Table 1). Panel c) plots the core-Sersic pn,,cS against the published Nukcr model's estimate of fib jvufc (Lauer et al. 2007b). Filled 
stars are the core galaxies from this study, while open circles are these same galaxies with R), and fn, from Lauer et al. 2007a. Triangles 
denote the seven Nuker 'core' galaxies which are reclassified here as Sersic galaxies. The solid lines (in panels a and b are Eqs. 9 and 10) 
are the least-squares fits to the fn,-L and fi^-cr core-Sersic data. Pearson correlation coefficients, r, and representative error bars, for our 
data, are shown at the bottom each panel. For the disk galaxies, the la uncertainty on the bulge magnitudes is ~ 0.75 mag. 



relation seen in Fig. [T7b . 

Fig. [Ek plots the core-Sersic fib versus the Nuker fib. In 
over-estimating the core radii, where 7' w 1/2, the Nukcr 
model underestimates the associated surface brightness, 
by typically 1 mag arcsec -2 and up to 2 mag arcsec -2 
with respect to the core-Sersic model. In Figs. IT9a and 
119b we show relations involving fib (the F-band surface 
brightness at the core-Sersic break radius and the Nuker 
'cusp radius', from Table 2 and Lauer et al. 2007a re- 
spectively) with My and a. Akin to Fig. ITTk . the (core- 
Sersic)-identified core galaxies display a correlation be- 
tween the core-Sersic fib and My, such that the OLS 
bisector regression analysis gives 

fi b = (-1.02±0.10)(My + 22) + (16.00 ± 0.20). (9) 

The OLS bisector fit to the core-Sersic fib and a yields 

a 



relation 



l0£ 



Rt 
pc 



= (1.01±0.69) log 



M 



BH 



10 9 Mr. 



(2.07 ±0.33). 

(12) 

Given the two widely used black hole mass estimators, 
the Mbh-o relation (Ferrarese & Merritt 2000; Gebhardt 
et al. 2000) and the Mbh-L relation (Marconi & Hunt 
2003), we can construct the R^Mbh relation from the 
Rb-cr and Rb-L relations in section [7T] to further investi- 
gate the core size and black hole connection. Combining 
the Mbh-v relation from Graham et al. (2011, their final 
entry in Table 2, acquired using only elliptical galaxies) 
with the Rb~o relation (Eq. 5) we can derive the new 
Rb-MBH relation 



fib = (11.61 ±1.60) log 



200 km s" 



#6 

pc 



log — = (1.03± 0.20) log 



/ M, 



BH 



-(2.08 ±0.22), 



(14.75 ±0.24). 

(10) 



7.2. Core size versus black hole mass 

Although the relation between the break radius and 
the supermassive black hole (SMBH) mass is less funda- 
mental than the relation between the central mass deficit 
and the SMBH mass (e.g. Graham 2004; Ferrarese et al. 
2006, Merritt 2006), de Ruiter et al. (2005) and Lauer 
et al. (2007a) nonetheless argue for the existence of a 
good correlation between the former. Shown in Fig. [20] 
is the core-Sersic break radius plotted against the black 
hole mass Mbh for 8 core galaxies with direct black hole 
mass measurements (see Graham 2008b and Graham et 
al. 2011). The bisector fit to R b and M B h for these 8 
galaxies gives 



V io 9 Mq 

(13) 

which is in remarkable agreement with Eq. 12. Com- 
bining the Rb-L relation (Eq. 6) with the Mbh~L rela- 
tion for predominantly massive spheroids from Graham 
(2007, his equation 19) converted to the y-band using 
B - V = 0.9 (Fukugita, Shimasaku & Ichikawa 1995), 
yields 



|£) =( , 45± 0, 9)1 „ g (J^_ 



log 



M 



(0.63±1.73) log ( j+ (2.03 ±0.78). 

(11) 

However, while considering only 7 of the 8 core galaxies, 
after excluding the only disk galaxy (NGC 3607), the re- 
gression analysis of Rb and Mbh yields the more certain 



(2.03 ±0.16). 

(14) 

While the Mbh-L relation is likely to be non-(log- 
linear), Graham 2012b, the log-linear relation presented 
by Graham (2007) is dominated by massive spheroids 
and thus a good representation of the Mbh~L relation 
for the 'core' galaxies. 

It is worth noting that the scatter in the direct Rb- 
MBH relation established using the 7 elliptical core 
galaxies with direct black hole mass measurements (Eq. 
12) is large. Although we only have a limited number 
of galaxies with a measured black hole mass and Eq. 12 
is somewhat driven by the highest mass black hole 8 , Eq. 

8 Removing the highest mass black hole and excluding the disk 
galaxy gives a slope of 1.63±0.75. 
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Fig. 20. — The core-Sersic break radius Rb,cS (Table 2) plotted 
as a function of the black hole mass Mbh for 8 core galaxies with 
direct black hole mass measurements, obtained from Graham 2008b 
and Graham et al. 2011 and adjusted to distances given in Table 1. 
The solid line (Eq. 11) is a least-squares fit to all 8 core galaxies, 
while the dashed line (Eq. 12) is a least-squares fit to the 7 core 
elliptical galaxies (open crosses) after excluding the only lenticular 
galaxy NGC 3607 (open diamond) having a Pearson correlation 
coefficient r=0.75 as shown in the plot. A representative error bar 
is shown at the bottom of the panel (see text for details). 

12 is consistent (over-lapping error bars) with the two 
inferred relations (Eqs 13 and 14). Moreover, in contrast 
to the discussion in Lauer et al. (2007a, their Eqs 20 and 
21), there is a good consistency among the deduced rela- 
tions (Eqs. 13 and 14). The updated relations are not (i) 
contaminated by the inclusion of galaxies without cores 
nor (ii) based on galaxy rather than bulge luminosity for 
the disk galaxies (an issue discussed by Graham 2008a, 
his section 6), and (hi) use core-Sersic break radii. 

7.3. Central mass deficit 

Galaxy merging is believed to be a common occurrence, 
responsible for the morphologies of elliptical galaxies. 
Growing evidence, based on extensive numerical experi- 
ment has indicated that the merger remnant of collisions 
between nearly equal mass spiral galaxies resemble early- 
type galaxies (e.g. Toomre & Toomrc 1972; Hernquist 
1993; Somerville & Primack 1999; Cox et al. 2006; Naab 
et al. 2006; Naab & Ostrikcr 2009). 'Dry' mergers of 
elliptical galaxies may subsequently build 'core' galaxies 
(e.g. Faber et al. 1997; Khochfer & Burkert 2003). 

In 'core' galaxies, the inner most stars are thought to 
have been ejected by inwardly spiraling binary SMBHs in 
the course of such dry, i.e. dissipationless, galaxy merg- 
ers, producing the observed central luminosity deficit, 
Ldef, relative to the inward extrapolation of the outer 
Sersic profile. Multiplying this deficit by the appropriate 
stellar mass-to-light ratio gives the central mass deficit, 
McLef- More specifically, scattering of stars from the 
galactic nuclei in three-body interactions (since there is 
no significant amount of gas which may render dynami- 
cal friction) is thought to be the avenue through which 
the coalescing binary black holes will 'harden' and be 
delivered to the galaxy center (e.g. Begelman, Bland- 
ford & Rees 1980; Nakano & Makino 1999; Milosavl- 
jevic & Merritt 2001). Using N-body simulations, Mer- 
ritt (2006) showed that core formation is a cumulative 



process, where the total central mass deficit after N dry 
mergers is such that M<w ~ 0-5NMbh- Gualandris & 
Merritt (2011) further discuss the long-term evolution of 
black hole binaries affecting the central stellar distribu- 
tion, and scouring out cores having radii bigger than the 
influence of the binaries. 

Here we have a) identified the galaxies with scoured 
out 'cores' (Table [2]) and b) quantified how their sizes 
scale with the final black hole mass Mbh (section 7.2). 
Graham's (2004) estimation of the central stellar mass 
deficit (Mdef) through the employment of the core-Sersic 
model yielded Mdef ~ Mbh, a result later confirmed 
by Ferrarese et al. (2006). This mass deficit is in ac- 
cord with hierarchical galaxy formation models, where 
luminous galaxies are a by-product of about 1 to 2 ma- 
jor dissipationless merger. For reference, based on ob- 
servations of galaxy pairs, Xu et al. (2011) found that 
since z = 1, massive galaxies have experienced 0.4 to 1.2 
major- mergers (see also Bell et al. 2004). 

In Fig. [21] we present tentative central stellar mass 
deficits obtained using our core-Sersic model parameters 
and Eq. A19 from Trujillo et al. (2004) plotted against 
the observed or predicted (Graham et al. 2011) super- 
massive black hole mass for each core galaxy. While R^, 
Hb and 7 are well constrained from our core-Sersic fits, 
R e and n may be less well constrained (section [4Tj and 
as such we caution about over-interpreting the results 
in Fig. [2U In passing we again note that the break ra- 
dius and mass deficit for NGC 3607 may be too large, 
based on HST/NIC2 data analysed by Richings et al. 
(2011). Here we simply remark that the mass deficits 
scatter around 0.5 to 4 times the central supermassivc 
black hole mass. Mdef /Mbh ratios less than 0.5 im- 
ply that a minor merger event may have taken place (in 
the absence of loss cone refilling). Note that we follow 
Graham (2004) and assume a F-band stellar mass-to- 
light ratio of ~ 3.5 (Worthey 1994) to compute the mass 
deficits. For further reference, using N-body simulations, 
Kulkarni & Loeb (2011), measure mass deficits which are 
up to 5 times the mass of the central black hole. 

Lauer et al. (2007a) proposed an alternate method to 
quantify the central stellar mass deficit in terms of a core 
with zero transition region breaking to an outer power- 
law profile with negative, logarithmic slope (3. However, 
as revealed by Graham et al. (2003), the Nuker model /3 
is not a robust quantity but varies with the radial range 
of the light profile that one tries to model. The range 
of (mass deficit)-to- (black hole mass) values in Fig. [5T] 
is notably less than the values in Lauer et al. (2007a) 
which were as high as ~ 19 at Mbh = 10 9 M Q (for black 
hole masses predicted using their Mbh~c relation). Kor- 
mendy et al. (2009, their figure 42) also reported central 
mass deficits that were some ^5 to 20 times greater than 
the central black hole mass, raising further doubts over 
their method of analyzing the galaxy light profiles (see 
the model-independent analysis by Hopkins & Hernquist 
2010). 

Finally, using Nuker model parameters, Giiltckin et al. 
(2011) report a break radius of 0".93 for NGC 4382 and 
a central mass deficit of 5.9 x 10 8 M Q . Our analysis yields 
a break radius roughly three times smaller (0".32) and 
a central mass deficit which is also three times smaller 
(1.8 x 1O 8 M ). 

Having refined the core galaxy sample in this study, in 
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Fig. 21. — Tentative central mass deficit (M de f) versus black hole 
mass (Mbh) for 32 core galaxies. We used Graham (2008b) and 
Graham et al. (2011) for direct supermassive black hole (SMBH) 
mass measurements of 8 core galaxies (open crosses), while the M-a 
relation presented in Graham ct al. (2011) was used for estimating 
the SMBH masses of the remaining 24 core galaxies (filled stars). 
A Representative error bar is shown at the bottom of the panel. 

a follow-up paper we intend to acquire a greater radial 
range of the light-profiles for these galaxies, enabling a 
better estimate of the outer Sersic parameters and thus 
the central mass deficits. This will allow us to check for a 
positive Md e f JMbh correlation with host spheroid mass 
and Mbhi tentatively seen in Fig. [2T1 



8. ADDITIONAL NUCLEAR COMPONENTS 

Additional nuclear light is detected in 12 (~31%) of the 
39 galaxies. Table 3 presents a comparison of detections 
of additional nuclear components from different studies; 
as done by Cote et al. (2006, their Table 3). Although 
the majority of our sample are core galaxies (~82%), 5 
of the 12 nucleated galaxies are Sersic galaxies. We find 
good agreement with the work of Cote et al. (2006) in 
assigning the additional nuclear components with only 
two exceptions (NGC 4472 and NGC 4552). NGC 4552 
has a point-source AGN (Renzini et al. 1995; Carollo et 
al. 1997; Cappellari et al. 1999) with a radio flux ^103 
mJy at 1.4 GHz (Condon et al. 1998). As for the Sy2 
galaxy NGC 4472, at odds with Lauer et al. (2005) but in 
agreement with Cote et al. (2006), Ravindranath et al. 
(2001) also did not detect any additional nuclear com- 
ponent (Soldatenkov, Vikhlinin & Pavlinsky 2003). We 
do however note that this galaxy's apparent point-source 
(which has a radio flux ^752 mJy at 1.4 GHz: White <fc 
Becker 1992) is very faint and the presence of dust can 
lead to such small irregularities in profiles acquired after 
image deconvolution. 

Lauer et al. (2005) identified the nuclear light excess as 
a central upward deviation from the best Nuker model 
fit to the host galaxy surface brightness profile. Their 
identification is not always consistent with our analysis. 
The relative faintness and extended nature of the addi- 



tional nuclear components in NGC 1374 and NGC 4458 
are viable explanations as to why our detections are at 
odds with Lauer et al. (2005). We also do not detect ad- 
ditional nuclear light in our data for NGC 1399 nor NGC 
4406. Ravindranath et al. (2001) and Cote et al. (2006) 
also noted the absence of central light excess in NGC 
4406, while Gebhardt et al. (2007) have also reported 
that NGC 1399 lacks any second component. Lyuben- 
ova et al. (2008) do however provide tentative evidence 
for the existence of a nuclear star cluster (or swallowed 
globular cluster) in NGC 1399. Lastly, Rest et al. (2001), 
after adopting a conservative central nuclei assigning ap- 
proach reported the absence of additional nuclear com- 
ponents in NGC 4365 and NGC 4478, while all successive 
studies, including ours, identified central light excess in 
these two galaxies. 

Of the 5 Sersic galaxies with additional nuclear compo- 
nent, only one, NGC 7213 (Veron-Cetty & Vcron 1988), 
has a central light excess associated with nonthermal 
emission from an AGN. We further note that the re- 
maining 4 nucleated Sersic galaxies have a nuclear disk, 
non-AGN central nuclei or both. In contrast, we note 
the presence of AGN in (at least) 5 of the 7 nucleated 
core galaxies: NGC 741 (Condon et al. 2002); NGC 
4278 (Younes et al. 2010); NGC 4472 (Diehl & Statler 
et al. 2008); NGC 4552 (Carollo et al. 1997; Cappel- 
lari et al. 1999) and NGC 5419 (Capetti & Balmaverde 
2005). Such prevalence is in accord with the studies by 
Balmaverde & Capetti (2006) and Richings et al. (2011), 
see also Pellegrini (2010). The central light excess in the 
remaining 2 core galaxies is not AGN related: NGC 4365 
has a stellar cluster (Carollo et al. 1997) and NGC 6876 
has a double optical nucleus possibly from an inclined 
disk (Lauer et al. 2002). 

9. CONCLUSIONS 

We have re-modelled the major-axis, surface bright- 
ness profiles of 39 alleged 'core' galaxies from Lauer et 
al. (2005), using Sersic and core-Sersic models. We have 
additionally and simultaneously accounted for the point 
sources and additional nuclear components that were ex- 
cluded by the Nuker analysis. Consistent with earlier 
published works, we found that the Sersic and core-Sersic 
models yield a robust representation of the underlying 
light distributions of Sersic and core-Sersic galaxies, re- 
spectively, all the way to the resolution limit. The typical 
rms residual scatter is 0.02 mag arcscc -2 . 

The main results of this work are: 

1. We have identified 7 of the 39 'core' galaxies from 
Lauer et al. (2005) to be Sersic galaxies which do not 
have partially depleted cores relative to the inward ex- 
trapolation of their outer Sersic light profile. This situ- 
ation tends to arise in galaxies and bulges fainter than 
My « —21 mag. Such galaxies with spheroid Sersic in- 
dex n < 3 or velocity dispersion a < 183 km s _1 arc 
not likely to have partially depleted stellar cores. 

2. We provide physical parameters (i?(,,^&i7) for the 
cores of 32 'core' galaxies , derived using the core-Sersic 
model. 

3. Due to noise or real small scale structure, non- 
parametric core size estimations obtained by locating 
the maximum of the second logarithmic derivative of the 
(non-smoothed) light profile, i.e. the point of greatest 
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TABLE 3 

Comparison of the detection of nuclear components from different studies 





Rest et al. 


Ravindranath 


Laucr et al. 


Cote et al. Our result 


Galaxy 


Profile Type (2001) 


ct al. (2001) 


(2005) 


(2006) 



NGC 0741 


c-S 






Yes 




Yes 


NGC 1374 


S 






No 




Yes 


NGC 1399 


c-S 






Yes 




No 


NGC 4278 


c-S 




Yes 


Yes 




Yes 


NGC 4365 


c-S 


No 




Yes 


Yes 


Yes 


NGC 4406 


c-S 




No 


Yes 


No 


No 


NGC 4458 


S 






No 


Yes 


Yes 


NGC 4472 


c-S 




No 


Yes 


No 


Yes 


NGC 4478 


S 


No 




Yes 


Possibly 


Yes 


NGC 4486B 


S 






No 


possibly 


Possibly 


NGC 4552 


c-S 






Yes 


No 


Yes 


NGC 5419 


c-S 






Yes 




Yes 


NGC 6876 


c-S 






No 




Possibly 


NGC 7213 


S 






Yes 




Yes 



curvature, appear to be unreliable (Fig. 

4. As with the Nuker model, the break radius of 
the core-Sersic model is shown to coincide with the ra- 
dius where it has a maximum in the second logarithmic 
derivative (Fig. 9, Right). 

5. For the first time, the radius where the negative 
logarithmic slope of the light profile 7' = 1/2, considered 
to be a suitable estimator for the size of the core, is 
shown to be consistent with the core-Sersic model break 
radius (Fig. ITTj). It should, however, be noted that even 
galaxies without depicted cores will have a radius where 
7' equals 1/2. Therefore, this measurement cannot be 
used to identify 'true' depleted-core radii. 

6. We have compared the core-Sersic break radii with 
the Nuker break radii. In line with previous works, we 
found that the Nuker break radii are larger than the 
core-Sersic break radii and also i? 7 < = i/2: on average, the 
Nuker break radii are ~2 times bigger than the core- 
Sersic break radii. Furthermore, the surface brightnesses 
(fib) at the Nuker model's break radii are up to 2 mag 
arcsec -2 fainter than the surface brightness at the core- 
Sersic model break radii. 

7. We have updated various structural parameter re- 
lations after excluding galaxies which do not have 'real' 
cores, and using core-Sersic parameters. We have also 
used the bulge magnitude instead of the galaxy magni- 
tude for the disk galaxies. We provide updated Rb-L, 
Rb-cr, Rb-^b, Hb~L and /i&-<7 relations in section 7.1. 
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Fig. 22. — Major-axis surface brightness profiles for the galaxies listed in Table 1. All profiles were obtained with the F555W (~ V-band) 
filter, except for NFC 1374 and NGC 1399 which were obtained with the F606W (~ ij-band) filter. The dashed curves show the Sersic 
component of the core-Sersic fits to the data, while additional nuclear sources were fit with either a Gaussian (dash-dot-dot-dot curve) or 
an exponential function (dotted curve). The solid curves show the complete fit to the profiles, with the rms residuals, A, about each fit 
given in the lower panels. Data points excluded from the fits are shown by the open circles. Brief detail on NGC 1374; NGC 3607; NGC 
3640 is provided in Appendix B, see also the text. 
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Fig. 22. — continued, see Appendix B text for further details. 
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Fig. 22. — continued, see Appendix B text for further details. 
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Fig. 23. — Similar to Fig. 1221 but for galaxies with various morphological peculiarities. See Appendix B text for further details. 
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12. APPENDIX B 

12.1. Notes on selected individual galaxies 

NGC 584: This galaxy has a depleted core with an inner negative logarithmic slope 7 = 0.46. Thus, it has 0.3 < 7 < 
0.5 as opposed to the core (7 < 0.3) / power-law (7 > 0.5) dichotomy (see also Glass et al. 2011). 
NGC 1374: This galaxy hosts an elongated inner disk (van der Marel & Franx 1993), which is well fit with a Sersic 
model plus an inner exponential function for representing the disk. 

NGC 3607: There is vivid evidence for the presence of an inner dust ring (Lauer et al. 2005). 

NGC 3640: The galaxy is known for its morphological peculiarity including several sharp and diffuse features (e.g. 
shells, ripples), which probably signal an ongoing or recent merger (Michard & Prugniel 2004), which makes the 
underlying host galaxy light difficult to model since our model is not designed to accommodate such peculiarities. An 
ongoing or recent merger could be the reason for having an unusually small core for a galaxy of this luminosity. 
NGC 3706: This galaxy exhibits an obvious inner stellar ring, shifting the peak of the surface brightness to R « 0".14 
from the photometric center (Lauer et al 2005). The light profile is better modelled with careful omission of the 
additional ring of starlight from 0".ll - 0".4 (see also Capetti & Balmaverde 2005). 

NGC 4073: A cD galaxy with a possible history of cannibalism and having similar distinct features as NGC 3706. 
Lauer et al. (2005) identified an asymmetric central ring of stars noticeable on the brightness profile of the galaxy. 
Data points from 0".l - 0".4 are excluded from the fit. 

NGC 4365: This object is a giant elliptical in the Virgo cluster with a kinematically distinct core (Forbes et al. 
1994), that has a slightly elongated inner nuclear excess (Carollo et al. 1997; Cote et al. 2006). The surface brightness 
profile is well fitted with the corc-Scrsic model plus a Gaussian function for the inner nucleus. 

NGC 4458: The isophotal contour analysis presented in Trujillo et al's. (2004) revealed a prominent nuclear disk, 

see also Morelli et al. (2004). Analysis of the brightness profile indicates the presence of an extended point source 

(Ferrarese et al. 2006) and an inner disk. A good match to the .ffST-observed profile is obtained fitting a Sersic model 

plus an inner exponential (n = 1) disk component and a Gaussian (n = 0.5) for the point source. 

NGC 4478: Like NGC 4458, there evidence for the presence of a nuclear disk in this galaxy (e.g. Trujillo et al. 

2004; Morelli et al. 2004). The only noticeable difference with the light profile of NGC 4458 is that this galaxy has a 

relatively elongated disk, and a compact point source (Carollo et al. 1997; Ferrarese et al. 2006). 

NGC 4486B: Based on the deconvolved HST /WFPC2 I- and V^-band images, Lauer et al. (1996) showed the existence 

of a double optical nucleus in this galaxy. However, the double nuclei are not obvious from most of the archival optical 
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HST images, as noted by Ferrarese et al. (2006), and it may be an artifact of the deconvolution routine. Nonetheless, 
the plateau in the inner light profile is not due to a relative deficit of stars. Tidal truncation from the interaction with 
the close companion M87 is apparent from the fit to the surface brightness profile, particularly in the outer part of the 
profile (R > 7" w 0.7 kpc). 

NGC 4552: This galaxy hosts a compact point source which is detectable in the optical image (Renzini et al. 1995; 
Carollo et al. 1997; Cappellari et al. 1999; Lauer et al. 2005). 

NGC 5419: There seem to be two compact nuclear point sources, the brighter one at the photometric center, visible 
in the optical images (e.g. Capetti & Balmaverde 2005; Lauer et al. 2005). A "dip", in the region 0".l to 0".2, is 
detected in the light profile of the galaxy. Thus, two data points from 0".l < R < 0".2 are excluded from the fit. 
NGC 6876: The dominant elliptical in the Pavo group shows past or ongoing interaction with the large spiral NGC 
6872 (Machacek et al. 2005). Furthermore, the archival I- and V-band images of this galaxy indicate the presence of 
a double optical nucleus (Fig. [4j, possibly the semi-digested nuclei of lesser galaxies or the ends of an inclined ring 
(Lauer et al. 2002). Like NGC 4486B, a core-Sersic model can fit the surface brightness profile (sec Fig. l2"5|) . Although, 
as the nuclei appear only 0".16, and equidistant, from the center, they may not explain the core's structure and break 
radius of 0.45". 

NGC 7213: This galaxy is the only spiral (Sa) in our sample. It hosts a bright Seyfert nucleus (Veron-Cetty & 
Veron 1988) and an inner disk. Hamced et al. (2001), see also Grosb0l et al. (2004), argue, based on the HI (neutral 
hydrogen) map, that this a highly disturbed system that may have experienced a past merger. In addition, using the 
broadband HST images, Deo, Crenshaw & Kraemer (2006) noted nuclear dust features in this galaxy. 



